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ABSTRACT 

We present forecasts for constraints on cosmological models which can be obtained by forth- 
coming radio continuum surveys: the wide surveys with the LOw Frequency ARray (LOFAR), 
Australian Square Kilometre Array Pathfinder (ASKAP) and the Westerbork Observations of 
the Deep APERTIF Northern sky (WODAN). We use simulated catalogues appropriate to the 
planned surveys to predict measurements obtained with the source auto-correlation, the cross- 
correlation between radio sources and CMB maps (the Integrated Sachs-Wolfe effect), the 
cross-correlation of radio sources with foreground objects due to cosmic magnification, and 
a joint analysis together with the CMB power spectrum and supernovae. We show that near 
future radio surveys will bring complementary measurements to other experiments, probing 
different cosmological volumes, and having different systematics. Our results show that the 
unprecedented sky coverage of these surveys combined should provide the most significant 
measurement yet of the Integrated Sachs-Wolfe effect. In addition, we show that using the 
ISW effect will significantly tighten constraints on modified gravity parameters, while the 
best measurements of dark energy models will come from galaxy auto-correlation function 
analyses. Using the combination of EMU and WODAN to provide a full sky survey, it will 
be possible to measure the dark energy parameters with an uncertainty of {cr(wo) = 0.05, 
c(wa) = 0.12} and the modified gravity parameters {cr(77o) = 0.10, (T(/io) = 0.05}, assum- 
ing Planck CMBH-SN(current data) priors. Finally, we show that radio surveys would detect a 
primordial non-Gaussianity of /nl = 8 at 1-(t and we briefly discuss other promising probes. 

Key words: large-scale structure of the universe — cosmological parameters — cosmology: 
observations — radio continuum: galaxies. 



1 INTRODUCTION 

Radio surveys for cosmology are entering a new phase with the 
construction of the LOw Frequency ARray (LOFAR, Rottgering 
2003), the Australian Square Kilometre Array Pathfinder (ASKAP, 
Johnston et al. 2008) and APERTIF, the new Phased Array Feed 
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receiver system for the Westerbork Synthesis Radio Telescope 
(WSRT, Oosterloo et al. 2010). In each case, the increased sensi- 
tivity available, together with a very wide sky coverage, will allow 
certain cosmological statistics to be measured with substantial ac- 
curacy. Several studies in the past have concentrated on the cosmo- 
logical constraints that can be determined from the large redshift 
surveys using the Hi 21-cm emission line (e.g. Abdalla & Rawl- 
ings 2005; Abdalla et al. 2010). However, little attention has been 
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paid to the information that can be gleaned from the large radio 
continuum surveys which will, in many respects, be much easier 
to interpret than the Hi surveys and allow us to push out to much 
higher redshifts. In this paper we consider three experiments us- 
ing the deep continuum observations, the auto-correlation of radio 
sources, the cross-correlation of radio sources with the Cosmic Mi- 
crowave Background (the late Integrated Sachs- Wolfe effect), and 
cross-correlation of radio sources with foreground objects (magni- 
fication bias). The level of the accuracy of these measurements, and 
the relative significance of the various potential probes, are the key 
issues which we wish to address in this paper. 

One of the goals of these measurements will be to measure the 
cosmological parameters of particular current interest. Among the 
biggest challenges in cosmology is to understand if the standard A 
Cold Dark Matter model, and its General Relativity context, is cor- 
rect, or if we need a different cosmological model and/or gravita- 
tional theory, with the related important implications for fundamen- 
tal physics. We will therefore present forecasts of the constraints 
on cosmological models and gravitational parameters that will be 
possible to obtain with the LOFAR, ASKAP and WSRT radio tele- 
scopes, in isolation and together. 

There are many major optical and near infra-red galaxy sur- 
veys (e.g. the Baryon Oscillation Spectroscopic Survey (BOSS, 
Eisenstein et al. 2011), BigBOSS (Schlegel et al. 2011), the Dark 
Energy Survey^ (DES), the Panoramic Survey Telescope And 
Rapid Response System^ (Pan-STARRS), Euclid (Laureijs 2009), 
the Large Synoptic Survey Telescope"^ (LSST)) which aim to im- 
prove the precision of cosmological parameter measurements dur- 
ing this decade. One of the goals of this paper is to discover whether 
there are also significant and complementary opportunities for im- 
provement of cosmological constraints by forthcoming radio con- 
tinuum surveys. These surveys have a niche because of their large 
sky coverage, high median redshift and number of objects ob- 
served. 

The paper is organised as follows. In Section 2 we will de- 
scribe the next generation of radio surveys, given by LOFAR, EMU 
and WODAN. In Section 3, we will discuss the predictions for 
source densities and bias as a function of redshift for each survey 
and for different source populations. In Section 4 we will present 
the cosmological probes we will use, and in Section 5 we show 
our predicted cosmological measurements. In Section 6 we will 
describe the methodology used to predict the resulting constraints 
on dark energy and modified gravity models, and in Section 7 we 
present our results. In Section 8 we present our conclusions and 
summarise why LOFAR, EMU and WODAN will be important for 
cosmology. 



2 FORTHCOMING RADIO SURVEYS 

In this section we introduce the three large radio surveys which 
we will focus on in this paper: LOFAR, EMU and WODAN. We 
do not consider the surveys to be conducted with the South African 
SKA Precursor Telescope (MeerKAT, Jonas 2009) as the parameter 
space probed by MeerKAT is towards much deeper and narrower 
surveys which are more adept to studying galaxy formation and 
evolution. In each case we will discuss the properties of the surveys. 
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^ http://www.lsst.org/lsst 



and their expected timescales for observation. A summary of the 
survey properties is shown in Table 1 . 

2.1 LOFAR 

LOFAR (the LOw Frequency ARray for radio astronomy, 
Rottgering 2003) is a multi-national telescope that has stations 
spanning Europe. The core of LOFAR is situated in the north-east 
of the Netherlands, with stations on longer baselines both within 
the Netherlands and across to Germany, UK, France and Sweden. 
Other stations may also be added throughout the rest of Europe in 
the coming years. 

Each LOFAR station operates at two broad frequency ranges, 
the high band operating at 120 < < 240 MHz and the low band 
which operates at 10 < < 80 MHz. The bulk of the early op- 
erations of LOFAR will be dedicated to a number of Key Sci- 
ence Projects (KSPs): Solar Physics and Space Weather, Transients, 
Cosmic Magnetism, the Epoch of Reionization, Cosmic Rays and 
Continuum Surveys. It is the last of these which is pertinent to the 
aims of this paper. 

The Continuum Surveys KSP will explore the bulk of the 
northern sky at low-radio frequencies. Low-frequency radio ob- 
servations are ideally suited for carrying out sensitive surveys of 
the extragalactic sky, firstly because the low-frequency ensures a 
large instantaneous field of view, allowing an increased survey 
speed compared to similar telescopes operating at higher frequen- 
cies. Second, the bulk of radio emission detected from extragalac- 
tic sources is due to synchrotron radiation and therefore increases 
towards the lower frequencies, although at the very lowest frequen- 
cies one might expect a turnover to occur due to synchrotron self 
absorption. 

The LOFAR continuum surveys (Rottgering 2010) follow the 
usual strategy of a "wedding cake" tiered survey. The design of 
these has focused on addressing the original key science topics 
within the continuum surveys, namely tracing the formation of 
massive galaxies, clusters and black holes using high-redshift ra- 
dio sources, measuring the star-formation history of the Universe 
through radio emission and tracing intracluster magnetic fields us- 
ing diffuse radio emission. However, as we demonstrate, these sur- 
veys will also provide key data which can be used to constrain the 
cosmology and gravitational physics in our Universe. 

For the purposes of this paper we concentrate solely on the 
120 MHz surveys from LOFAR, as these are the most sensitive for 
our science, i.e. wide-field and highly sensitive. The tiers of the 
LOFAR survey are as follows: the large-area, "Tier-l" survey will 
survey the whole of the northern sky down to an expected rms flux- 
density at 120 MHz of 5120 mhz = 0.1 mjy. 

The LOFAR Tier-2 survey will survey to deeper levels over a 
smaller area. The baseline strategy is to survey around 550 square 
degrees at 120 MHz to an rms flux-density of 5*120 mhz = 25/iJy. 

The Tier-3 survey is not considered in this paper due to the rel- 
atively small area it will survey (~ 70 square degrees) at 150 MHz 
to ~ 6/^Jy rms. 

We also consider what results could be achieved with the LO- 
FAR commissioning survey. Although this survey is still being fully 
defined, we take a shallow survey covering the whole northern 
hemisphere at 150 MHz and use a lOcr limit of 7 mJy. This al- 
lows us to examine what can be achieved with a very conservative 
survey. 

If we assume only the stations situated in the Netherlands are 
used in carrying out the large area surveys, then the resolution at 
120 MHz will be ~ 6 arcsec and around ~ 5 arcsec at 150 MHz, 
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Survey 


Area 


Frequency 




Mean z 


Median z 


LOFAR MS^ 


2tx 


150 MHz 


1.0 X 10" 


1.6 


1.3 


LOFAR Tierl 


2-K 


120 MHz 


6.5 X 10** 


1.8 


1.1 


EMU 


Stt 


1400 MHz 


2.2 X 10'' 


1.7 


1.1 


WODAN 


Itt 


1400 MHz 


7.3 X 10" 


1.7 


1.1 



Table 1. Parameters of the surveys considered. We use the lOcr flux-density 
limit for each survey. Total number of radio sources Ng^i, mean and median 
redshifts calculated using our number density models in Section 3. 



3 SOURCE POPULATION MODELS 

In this section we describe our models for source populations for 
LOFAR, EMU and WODAN surveys; in particular the number den- 
sity of different source populations as a function of redshift, and the 
bias of different source populations as a function of redshift. These 
are required in order to make predictions for cosmological probes 
such as the auto correlation function, the ISW effect and magnifi- 
cation bias. 



i.e. very similar to the Faint Images of the Radio Sky at 21 centime- 
tres (FIRST) survey (Becker et al. 1995). 



2.2 EMU 

EMU (Evolutionary Map of the Universe, Norris et al. 2011) is 
an all-sky continuum survey planned for the new Australian SKA 
Pathfinder (ASKAP) (Johnston et al. 2008) telescope under con- 
struction on the Australian candidate SKA site in Western Aus- 
tralia. EMU is one of the two key projects (the other is the WAL- 
LABY all-sky HI survey) which are primarily driving the ASKAP 
design. At its completion, expected to be in late 2012, ASKAP will 
consist of 36 12-metre antennas spread over a region 6 km in diam- 
eter. Although the collecting area is no larger than many existing 
radio telescopes, the phased- array feed at the focus of each an- 
tenna provides about 100 dual-polarisation pixels, giving ASKAP 
a thirty square degree of instantaneous field of view. This enables 
it to survey the sky some thirty times faster than existing radio tele- 
scopes at similar frequencies. 

The primary goal of EMU is to make a deep (10 /iJy rms) 
radio continuum survey of the entire Southern Sky, extending as far 
North as +30 deg. EMU will cover the same area (75% of the sky) 
as NVSS (Condon et al. 1998), but will be 45 times more sensitive, 
and will have an angular resolution (10 arcsec) five times better. 
It will also have higher sensitivity to extended structures. EMU is 
expected to begin in late 2012 and it will generate a catalogue of 
radio sources 38 times greater than NVSS; all radio data from the 
EMU survey will be placed in the public domain as soon as the data 
quality has been checked. 



2.3 WODAN 

WODAN (the Westerbork Observations of the Deep Apertif North- 
ern sky survey) is planned to chart the entire northern sky above 
Dec > 30° down to a proposed rms flux density at 1.4 GHz of 
51.4 GHz = 10/iJy/beam (Rottgering et al. 2011). It will be able to 
do this because of the new phased array feeds (APERTIF) being put 
on the Westerbork Synthesis Radio Telescope (WSRT, Oosterloo et 
al. 2010). The phased array feeds will open up the field-of-view of 
the WSRT to around 8 square degrees allowing very high survey 
speeds. Such a survey is extremely complementary to the proposed 
LOFAR Tier- 1 survey and will allow source spectral indices to be 
measured down to very faint levels. Although APERTIF increases 
the field-of-view of the WSRT considerably it will remain a rela- 
tively low-resolution survey instrument, with the resolution limited 
to the distribution of the WSRT antennae; as such the resolution 
will be around ~ 15 arcsec. However, this resolution is generally 
not a problem for the experiments we discuss in this paper. The cur- 
rent schedule for the commencement APERTIF surveys is 2013. 



3.1 Number densities 

We use empirical simulations to predict the number density of ra- 
dio sources per redshift interval for the envisaged all-hemisphere 
LOFAR survey, the WODAN survey, and the ASKAP-EMU Stt 
steradian survey. The combination of these surveys will provide 
complete coverage of 47r steradians of the sky; however, the differ- 
ent observing frequencies and depths means that they will produce 
distinct redshift distributions, which need to be understood in order 
to use the combination for cosmological constraints. Throughout 
this paper we assume that no redshift information is available for 
individual radio sources. 

We use the simulations of Wilman et al. (2008, 2010), de- 
veloped for predictions for the Square Kilometre An^ay continuum 
survey. These simulations provide specific prescriptions for the red- 
shift evolution of the various populations which dominate the radio 
source counts: powerful active galactic nuclei at bright fluxes, down 
to the less luminous radio-quiet AGN, starburst and star-forming 
galaxies. The simulations cover five different radio frequencies - 
150, 610, 1400, 4860 and 18000 MHz. We use flie update of the 
simulated catalogue (Wilman et al. 2010), which has been adjusted 
to incorporate results from mid- and far-infrared data to provide a 
better estimate of the starburst and star-forming galaxy populations. 

The N(z) from these simulations should in principle be modi- 
fied by Redshift-Space Distortions (Rassat 2009) and magnification 
bias (Loverde et al. 2007), which are not included. However these 
corrections are small and will not affect our results. 

Catalogues are generated from the S'^ database* correspond- 
ing to the radio flux-density limits of the proposed LOFAR, EMU 
and WODAN surveys. As described in Section 2, we assume the 
depth of the LOFAR survey over the whole hemisphere to be uni- 
form across the sky down to a rms flux-density of 0.1 mjy as given 
in the LOFAR Surveys document Morganti et al. (2010). For the 
purposes of this paper we use the 151 MHz data from the S'^ data 
base and extrapolate to 120 MHz using the spectral index deter- 
mined between 610 MHz and 151 MHz to predict the number den- 
sity distribution. As the simulations include a spectral curvature 
term, this means that the spectral index between 610 and 151 MHz 
is generally flatter than the canonical a ~ 0.7 and thus leads to a 
decrease in the number of sources one would expect based on using 
an a ~ 0.7. We adopt these numbers as a conservative approach, 
however we note that if the radio spectra do not flatten significantly 
to low frequencies then the constraints from LOFAR Tier 1 will be 
similar to those of the EMU survey, albeit over 2 tt sr rather than 
3 TT sr. We then apply a cut to the simulated data and retain only 
those sources with an integrated flux density larger then 10 times 
the rms noise in the map. Note that for extended sources the nom- 
inal sensitivity to peak flux densities is then less than 10 a. How- 
ever, this definition ensures that the virtually all extended sources 

* http://s-cubed.physics.ox.ac.uk 
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will be detected. It is also conservative in the sense that it still has to 
be proven that all these new instruments can reach their theoretical 
thermal noise levels. 

For the EMU and WODAN surveys we again use the Wilman 
et al. (2008, 2010) simulations, this time at the 1.4 GHz frequency. 
EMU will survey ~ 75% of the sky down to an rms flux-density 
limit of 10/i Jy, while WODAN will survey 10,000 sq. deg. down 
to an rms flux-density limit of lO/i Jy; we again extract a catalogue 
from the S*^ database down to this limit and apply cuts at 10 u 
signal-to-noise. Note that Norris et al. (201 1) assume a less conser- 
vative 5 (T threshold and therefore obtain stronger constraints. 

In Figures 1 and 2 we show the resulting redshift distributions 
adopted for the different surveys; in Figure 1 we display the total 
number of radio sources for the LOFAR (MS^ and Tierl), EMU 
and WODAN surveys, while in Figure 2 we display the number of 
sources for the different source types (Star Forming Galaxies, Star- 
Bursts, Radio Quiet Quasars, Fanaroff-Riley 1 and Fanaroff-Riley2; 
see Wilman et al. (2008) for details of how each of these is defined) 
within the surveys. 

In addition to predictions for the four surveys considered, we 
will also consider a combination of EMU and WODAN, given that 
they will span a similar range of frequencies and depth (see Tab. 
1), in order to have a complete full sky catalogue covering both 
hemispheres. 

3.2 Galaxy Bias 

As we will often be using radio sources as a probe of large-scale 
structure, it is necessary to model how biased the sources are in 
relation to the underlying structures. On large scales we assume 
that the two-point correlation function can be written (Matarrese et 
al. 1997, Moscardini et al. 1998) as: 



^(r, 2) = 6^(Mcff,2)CDM(r, 2) 



(1) 



where McB represents the effective mass of dark matter halos in 
which sources reside and ^da/ is the correlation function of dark 
matter. We derive a model of the bias using the peak-background 
split formalism (Cole & Kaiser 1989, Mo & White 1996), following 
the prescription of Sheth & Tormen (1999); in this context the mass 
function of halos, altered from Press & Schechter (1974), is given 
by: 



n{z,M) 



TV 8-kG MD{z)aM 



(2) 



1 + 



■ exp 



fD{z)aMy''' 




din am 


V V^Sc J 




din A/ 



where aJ^J is the mass variance on scale M, 5c is the critical over- 
density for the spherical collapse, D[z) the linear growth factor of 
density fluctuations, and q and p are parameters to be fitted with 
simulations (Sheth & Tormen 1999); we can describe the bias by: 



h{M,z) + Y 

Oc 

5c 



q5l 



D2{z) 



+ 



(3) 



l + {V^Sc/[D{z)aM]r^ 



For the purposes of this paper we use the bias in the 5*^ simulation 
for each galaxy population, which is computed using the formalism 
of Eq. (2), (3) separately for each galaxy population, where each 
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Figure 3. Bias as a function of redshift for the different source types, as 
calculated for our simulated catalogues in accordance with Wilman et al. 
(2008). 



population is assigned a dark matter halo mass. This dark matter 
halo mass is chosen to reflect the large-scale clustering found by 
observations. Note that for most of Fig. 3, there are simply no ob- 
servational measurements available at present, so large uncertain- 
ties in bias remain. 

The simulation provides us with a source catalogue with 
the sources identified by type, i.e. starburst, FRII-type radio galaxy 
etc. Each of these has a different prescription for the bias, as de- 
scribed in Wilman et al. (2008). With this framework, one finds 
that the increasing bias 6(2) with redshift would lead to exces- 
sively strong clustering at high redshift, therefore the bias for each 
population is held constant above a certain cut-off redshift, as de- 
scribed by Wilman et al. (2008). The resulting redshift dependence 
of the bias we use for the different source types is shown in Fig. 3. 
While this bias evolution is indicative of that expected, the exact be- 
haviour is not yet well known; to allow for this uncertainty, we will 
marginalise over the overall bias amplitude, and discuss remaining 
uncertainties in Section 8.1. 



4 COSMOLOGICAL PROBES 

In this section, we will describe several cosmological probes which 
one can measure with the forthcoming radio surveys, in combi- 
nation with surveys at other wavelengths. Here we describe the 
necessary framework for calculating the accuracy with which we 
can measure these probes; in the next section we will describe the 
signal-to-noise of measurements with the specific planned surveys. 



4.1 The Auto-Correlation Function 

The first cosmological probe we can examine with the radio surveys 
is the two-point correlation function of source positions, which is 
a measure of the degree of clustering in either the spatial, ^(r), or 
the angular, distribution of sources. For the current radio sur- 
veys, where individual redshifts will be unknown, there will be lit- 
tle radial information, therefore it will be appropriate to only study 
the angular two-point correlation function, w{Q\ which is defined 
as the excess probability of finding a radio source at an angular 
distance Q from another given radio source (Peebles 1980): 



5P = n\V^w{d)\5Q,, 



(4) 
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where 5P is the probabiUty, n is the mean surface density and 
a surface area element. 

The angular two-point correlation function of a given sample 
of objects can be computed using one of the many estimators that 
have been proposed (e.g. Hamilton 1993, Landy & Szalay 1993). 

Its Fourier transform, the angular power spectrum, can be cal- 
culated from the underlying 3D matter power spectrum using: 

cr = K,„a?;> = 4^/ f A^(fc)[ly/(fc)l^ (S) 

where is the radio source distribution window function de- 
scribed below in Eq. 7, and A^(fc) is the logarithmic matter power 
spectrum today, and aim are the spherical harmonics coefficients, 
assumed to be standard gaussian random variables. 

4.1.1 The radio source window function 

Radio source counts are a biased tracer of the underlying matter 
distribution, and thus the projected number density of radio sources 
per steradian is related to the matter distribution via: 

dN 

n(z,h)dzdQ = [1 + b(z)6(z,n)]dzdQ. (6) 

dz 

The window function can then be written as (see e.g. Raccanelli et 
al. 2008; Giannantonio et al. 2008a): 

W^ik) = f ^b{z)D{z)je[ckij{z)]dz, (7) 



where {dN/dz)dz is the mean number of sources per steradian 
with redshift z within dz, brighter than the flux limit, b{z) is the 
bias factor relating the source overdensity to the mass overdensity, 
assumed to be scale-independent, D{z) is the linear growth factor 
of mass fluctuations, ji{x) is the spherical Bessel function of order 
£, and rj(z) is the conformal look-back time. 

4.1.2 Non Gaussian clustering 

The amplitude and shape of clustering on large scales, described 
by the ACF, can provide important cosmological information. For 
example, a unique way to test aspects of inflationary theories is 
given by measuring the statistics of the initial conditions of cos- 
mological perturbations. An important goal for forthcoming cos- 
mological experiments is to test whether initial conditions of the 
probability distribution function of cosmological perturbations de- 
viate from gaussianity; this can be done using the CMB (Bartolo et 
al. 2004, Komatsu 2010 and references therein) or the large-scale 
structure of the Universe (Matamse et al. 2000, Dalai et al. 2008, 
Slosar et al. 2008, Desjacques & Seljak 2010, Xia et al. 2010). De- 
viations from Gaussian initial conditions can be parametrized by 
the dimensionless parameter /nl: 

^NG = </> + /nl (<?!>' - ((/>')) , (8) 

where $ denotes Bardeen's gauge-invariant potential, which, on 
sub-Hubble scales reduces to the usual Newtonian peculiar gravi- 
tational potential. Here (f) is a Gaussian random field, and the sec- 
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Figure 2. Redshift distributions of LOFAR, EMU and WODAN surveys, for different source types: Star-Forming Galaxies, Star Burst, Radio Quiet Quasars, 
FRl and FR2 sources; on vertical axes are the number of sources per bin of width Az = 0.3. 



ond term, when /nl is not zero, gives the deviation from gaussian- 
ity; in this paper we refer to the so-called "local type" /nl and 
we use the LSS convention (as opposed to the CMB one, where 
/n^ ~ 1-3 fS^"", (Xia et al. 2010)). 

One method for constraining non-Gaussianity from large- 
scale structure surveys exploits the fact that a positive /nl corre- 
sponds to positive skewness of the density probability distribution, 
and hence an increased number of massive objects (Matarrese et al. 
2000, Dalai et al. 2008, Desjacques & Seljak 2010). 

In particular, a non-zero /nl in Eq. (8) introduces a scale- 
dependent modification of the large-scale halo bias, so that the dif- 
ference from the usual Gaussian bias, is: 

AK.,.)^[.o(.)-l]/.L4. ^,g-g(^^ , (9) 

where boiz) is the usual bias calculated assuming gaussian ini- 
tial conditions, assumed to be scale-independent, D{z) is the linear 
growth factor and Sec is the critical value of the matter overdensity 
for ellipsoidal collapse. Sec = Sc^/q". 

4.2 The Integrated Sachs-Wolfe Effect 

In addition to making an auto-correlation of source positions, it is 
possible to cross-correlate the radio source distribution with CMB 
temperature maps, in order to detect the so-called Integrated Sachs- 
Wolfe (ISW) effect (Sachs & Wolfe 1967). Travelling from the last 
scattering surface to us, CMB photons pass through gravitational 



potential wells of intervening matter. In an Einstein-de Sitter uni- 
verse, the blueshift of a photon falling into a well is cancelled by the 
redshift as it climbs out. However, in a universe with a dark energy 
component or modification to General Relativity, the local gravita- 
tional potential $ varies with time, so potential wells are stretched 
while photons are traversing the well; this leads to a net blue-shift 
of the photons, and equivalently to a net change in photon tempera- 
ture, which accumulates along the photon path, and is proportional 
to the time variation of the gravitational potential. 

The integrated Sachs-Wolfe effect only contributes to the low 
I multipoles of the CMB fluctuations, and is smaller than the pri- 
mary CMB anisotropics even at those I. Thus, to make the ef- 
fect detectable, we have to cross-correlate CMB maps with tracers 
of large scale structure (Crittenden & Turok 1996) such as radio 
sources, since the source density traces the potential wells. If the 
evolution of potentials is modified by dark energy we should ob- 
serve a correlation between CMB temperature anisotropics and the 
source distribution; for this reason, ISW measurements will provide 
a signature for dark energy or modified gravity. The WMAP data 
have been cross correlated with a variety of radio, IR, optical, and 
X-ray surveys (e.g. Giannantonio et al. 2006, 2008a, Pietrobon et 
al. 2006, Raccanelli et al. 2008; see Dupe et al. 2010 for a review of 
recent results and more references) to look for evidence of a decay 
of the gravitational potential due to the influence of dark energy. 

We can write the cross-correlation power spectrum between 
the surface density fluctuations of radio sources and CMB temper- 
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ature fluctuations as: 

Cf = =4^/ '^A\k)W!{k)Wl{k), (10) 

where Wl and Wj are the radio source and CMB window func- 
tions, respectively, and (fc) is the logarithmic matter power spec- 
trum today. 

The cross-correlation function as a function of the angular sep- 
aration 9 is then obtained as: 



(11) 



where Li are the Legendre polynomials of order I (Legendre 1785). 



4.2.1 The ISW window function 

In the Newtonian gauge, scalar metric perturbations are specified 
by the gauge-invariant potentials ^ and $: 



ds 



-a^(r)[(l + 2*)dr^ 



2^)dx 



(12) 



The temperature anisotropies due to the ISW effect are expressed 
by an integral over the conformal lookback time from today (77 = 
0) to the CMB decoupling surface rjdec'- 

eisw^ — = -- ($ + *)d77, (13) 

where r is the conformal time, the dot represents a conformal time 
derivative and the integral is calculated along the line of sight of 
the photon. In the absence of anisotropic stress, the momentum 
constraint in GR fixes "I> = —5', so the ISW modification of the 
temperature of the CMB in GR becomes: 



ST 2 r"" , 
Qisw = — = I —dV- 



(14) 



10 dri 

The local gravitational potential is related to the matter distri 
bution via the Poisson equation: 



(15) 



where the gradient is taken with respect to comoving coordinates; 
taking the Fourier transform we have: 



(16) 



where Hq is the Hubble constant, g{rj) = D{rj)/a{rj) is the linear 
growth suppression factor and S{k) is the mass overdensity field. 

Combining Eq. 15 and 16, the window function for the ISW 
effect can be written as: 

W^{k) = 3Qom (^^) J ^jl[ck7^{z)]dz, (17) 

where $(2:) is the Newtonian gravitational potential. 



4.3 Magnification bias 

Light rays are deflected by large scale structures along the line of 
sight, which therefore systematically introduce distortions in the 
observed images of distant sources; this is the phenomenon of grav- 
itational lensing. The sources behind a lens are magnified in size, 
while surface brightness is conserved; this leads to an increase in 
the total observed luminosity of a source. 

Observationally we can detect the effects of magnification by 



cross-correlating two galaxy surveys with disjoint redshift distribu- 
tions; in this paper, we consider the possibility of using an optical 
survey such as SDSS-II (Abazajian et al. 2009) or DES^ for our 
low redshift "lens" sample (Pan-STARRS will also be available on 
these timescales and could also be used), which will serve as the 
foreground which magnifies the background radio distribution; we 
will discuss this further in Section 5.3. 

This "cosmic magnification" effect was first detected by 
Scranton et al. (2005), who cross-corelated foreground SDSS 
galaxies with SDSS quasars. More recently, Hildebrandt et al. 
(2009) have detected the effect in samples of normal galaxies in 
the Canada-France-Hawaii-Telescope Legacy Survey, Wang et al. 
(2011) have detected the effect at longer wavelengths using Her- 
schel, while Menard et al. (2010) have built on the SDSS analysis 
by constraining galaxy-mass and galaxy-dust correlation functions. 

The effect can be described in detail as follows. At position 
if, we can relate the behaviour of unlensed sources with number 
density No{m)dm within a magnitude range [m, m + dm], to that 
of lensed sources with number density N{m, (p)dm. There are two 
competing effects in this relationship, namely the flux increase due 
to magnification of distant faint sources, which increases the num- 
ber density of observed images above a certain magnitude thresh- 
old; and counteracting this, the dilution of the number density due 
to the stretching of the solid angle by lensing. If the source fluxes 
have a distribution with a power law slope given by: 



a{m) = 2.5 



9[logjVo(m)] 
dm 



one can obtain (Bartelmann & Schneider 2001): 

N{m,(p)dm = ^°''"'~'^A''o(m)dm, 
where the magnification /i is: 



1(1-^)2 -lipi- 



ds) 



(19) 



(20) 



where the convergence k and the shear 7 are two further lensing 
distortions. In the weak lensing regime it is possible to Taylor ex- 
pand the last equality in Eq. (20) to obtain: 



^{(p) ~ 1 + 2k{0). 



(21) 



We therefore see that the magnification is closely related to the con- 
vergence K, which is related to the matter overdensity via a line-of- 
sight integral (Bartelmann & Schneider 2001): 



d'wW{w)fK{w 



5{fK{w)(p,w) 
a{w) ' 



(22) 



with (p being the angular position on the sky, wh the horizon dis- 
tance, 'w{z) the comoving radial distance, fK{w) the angular di- 
ameter comoving distance, a the scale factor, and a quantity W 
involving the redshift distribution and geometry: 



fK{w') 



(23) 



for which Zw{'w)dw is the source redshift distribution. 

Because of the magnification bias effect described, we can ob- 
tain cosmological constraints by cross-correlating foreground and 
background objects, and hence investigating how clustered lensed 
background sources appear to be around foreground sources, com- 
pared to a random distribution. The most common estimator of the 
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angular two point correlation (also adopted in this work) is given 
by: 



isL = [NsNLr^[Ns{v) - Ns][NL{(p + cj)) - Nl 



(24) 



where in our case 5* and L indexes denote background sources and 
foreground lenses and overbarred quantities correspond to averaged 
quantities. 

Large scale structures only slightly magnify or demagnify 
sources, so we can write: 



a — 1 /-I 



-S^l)°''' ~ 1 + {a-l)5ii, 
leading to an over/under-density in background sources: 

Ns{0)-Ns 



Ns 



l)5^l(^). 



(25) 



(26) 



Assuming that foreground sources have bias fe_L, the number density 
can be related to the underlying matter density contrast by: 



Nl 



= bL5{(p), 



(27) 



Then as a consequence (c.f. Bartelmann & Schneider 2001), 
^sl{'p) is related to the theoretical magnification density contrast 
2-point correlation function (<^) via: 



1)&l(<^)Cm4'?)> 



(28) 



with 



dw 



W{w)Gf{w) 
a{w)fj.{w) 



j kdkPs{k)Jo{kip), 



(29) 



where Ps{k) is the matter power spectrum, Gf{w) is the fore- 
ground redshift distribution and W{'w) is the source lensing effi- 
ciency distribution given in equation (23). 

It is only when a ^ 1 that we obtain magnification bias. We 
obtain a positive cross-correlation only when a > 1 and anticor- 
relation when a < 1. Equations (28) and (29) show how magni- 
fication bias observations allow us to measure information about 
the amplitude, shape and evolution of the matter power spectrum, 
together with information about the bias and geometrical factors in 
the expanding background. 

From (Eq. 29) we can also obtain the cosmic magnification 
power spectrum (Bartelmann & Schneider 2001): 

dk 



cr 



k 



■A\k)W^{k)Wt{k), 



(30) 



where has the same meaning as in Eq. (5), and W^{k) 

contains the prefactors and w integral from equations 28 and 29. 



5 PREDICTIONS FOR MEASUREMENTS WITH 
FORTHCOMING SURVEYS 

5.1 Autocorrelation predictions 

We computed the predicted auto-correlation source power spectra 
for LOFAR, EMU and WODAN using Eq. 5; the errors were as- 
sumed to follow: 



i2£ + l)/,ky 



(31) 



where /sky is the sky coverage of the survey and n is the mean 
number of sources per steradian. This assumes that systematics are 




Figure 5. Source power spectrum of EMU radio sources for different values 
of the non-Gaussianity parameter /nl: shaded regions are errors for the 
EMU survey as in Eq. 3 1 . 



sub-dominant, and there are no effect from the finite size of objects 
(i.e. we are not close to the confusion limit). 

In Fig. 4 we show for the combined source populations of 
the four different surveys considered. As shown, the errors on large 
and small scales are more pronounced since they are dominated 
by cosmic variance and shot noise, respectively; at intermediate 
scales, surveys with higher number density will provide the best 
measurements (in this case EMU). We will examine what can be 
learned cosmologically from these measurements in section 7. 

In Fig. 5 we plot the predicted source power spectrum of EMU 
radio sources for different values of the non-Gaussianity parameter 
/nl; the black solid line is the standard Gaussian prediction, the 
other lines being the prediction for non-Gaussian clustering, and 
the shaded area is l-a errors (per mode) as in Eq. 31. The presence 
of the non-Gaussian bias of Eq. 9 enhances the clustering at large 
scales, thus increasing the amplitude of the autocorrelation function 
at those scales. A analysis shows that EMU should be able to 
distinguish (at 1-cr level) a /nl of 8 from a purely Gaussian model; 
it is worth noting that the current limit on /nl from WMAP 7-year 
data is /^Jp = 42± 27 at 68% CL (Komatsu et al. 201 1*^), and any 
detection of /nl ^ 1 would rule out all single scalar field inflation 
models (Komatsu 2010). 

It is intriguing to note that the observed autocorrelation func- 
tion of the NVSS Survey has a shape that differs from the ACDM 
prediction at relatively large angular separation (Xia et al. 2010); 
the observed behaviour can be explained using a non-Gaussian cor- 
rection (Xia et al. 2010), a peculiar bias model (Raccanelli et al. 
2008), or some systematic errors not yet found in the NVSS Survey. 
The degeneracy between models of bias and non-Gaussian correc- 
tions can be broken because their effect has a different redshift and 
scale dependence (the non-Gaussian correction is important only at 
large scales, because of the 1/fc^ term in Eq. 9). 

It is also interesting to note that a similar excess power at large 
scales has been found in spectroscopic (Kazin et al. 2010, Samushia 
et al. 201 1) and photometric (Thomas, Abdalla & Lahav 2010) data 
sets, although Ross et al. (201 1) suggest that this is likely to be due 
to masking effects from stellar sources. With the forthcoming all- 
sky radio surveys, measuring the angular autocorrelation will be 



Note that in the CMB convention this is 32 ± 21 
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Figure 4. Source power spectra (Eq. 5) of the combined source populations (black solid lines) for the different surveys, with 1-cr errors (grey shaded regions), 
as in Eq. 32. 



an interesting check for this problem. On the other hand, the CMB 
shows a lack of correlation at angular scales > 60 degrees (Ben- 
nett et al. 201 1, Copi et al. 2010), discrepant with the concordance 
model of cosmology. The significance and origin of this is unclear. 
Certainly radio surveys that cover large fractions of the full sky will 
help to resolve this puzzle. 



5.2 Cross-correlation predictions 

As we have seen in section 4.2, the cross-correlation between the 
CMB and the LSS depends on various factors from both the win- 
dow functions in Eq. (10); it is influenced by the evolution of the 
gravitational potential (Eq. 17) and by the clustering and bias of 
structures (Eq. 7), and for this reason it has been used to test and 
constrain cosmological issues such as the evolution and clustering 
of structures (Raccanelli et al. 2008, Massardi et al. 2010, Schaefer 
et al. 2009), models of dark energy (Pogosian et al. 2005, Xia et al. 
2009) and alternative models for the gravitational potential, such as 
the DGP, Unified Dark Matter cosmologies and Brans-Dicke theo- 
ries (Giannantonio et al. 2008b, Bertacca et al. 2011, De Felice et 
al. in preparation). 

The detection of the ISW effect via the cross-correlation of the 
LSS with the CMB is cosmic-variance limited, as it affects only the 
largest angular scales; therefore, the best measurement possible is 
a complete full sky survey with negligible shot noise. Regarding 
the CMB, the data provided by WMAP is already precise enough 



at low i, and the improvement that Planck will provide does not 
substantially affect the ISW detection significance. 

In Figure 7 we show the predicted cross-correlations of the 
CMB with the combined radio source distributions. Solid lines are 
standard ACDM+GR model, shaded regions are errors, calculated 
via (see e.g. Cabre et al. 2007): 



{21 + l)/,ky 



(32) 



where /sky is the sky coverage of the survey. The shot noise should 
be negligible for the cross-correlation of the CMB with LOFAR 
Tierl, EMU and WODAN on the scales of interest, given the high 
number density per I mode. 

Fig. 8 shows the predicted cross-correlation function, with 
cosmic variance errors, for the WODAN survey (light grey area), 
the WODAN+EMU combination (dark grey area) and the mea- 
sured NVSS errors (error bars) as a comparison. Note the substan- 
tial improvement in ISW measurements provided by the all-sky sur- 
vey, compared to NVSS or WODAN alone. 

The enhanced clustering due to non-Gaussianity would also 
modify the cross-correlation of galaxies with the CMB (Xia et al. 
2010), through the modified bias of Eq. (9) in the galaxy window 
function (Eq. 7). The effect is more significant on the largest scales, 
as shown in Fig. 6, which presents the cross-correlation of the CMB 
with EMU radio sources for different values of /nl- As in the auto- 
correlation case, we performed a analysis to predict what level 
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Figure 6. Cross-power spectrum of EMU radio sources with tlie CMB 
(Eq. 10) for different values of the non-Gaussianity parameter /nl; shaded 
regions are en'ors for the EMU survey as in Eq. 32. 
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Figure 8. Cross-correlation of WODAN sources with the CMB. Black solid 
line is the ACDM prediction, black dashed line is the UDM prediction for 
= 10~^ (see text for details); shaded regions are en'ors, light grey for 
the WODAN survey, dark grey for the EMU-l- WODAN combination; eiTor 
bars are NVSS errors. 



of non-Gaussianity we should be able to detect; we used again the 
simulations of EMU data as its ISW detection should be the best 
of the surveys analysed (see Fig. 7), and we found that these data 
would detect a /nl of 11 at 1-cr level. 

As an initial example of cosmological constraints, in Fig. 8, 
the black dashed line is the predicted cross-correlation function for 
Unified Dark Matter scalar field cosmologies, where dark matter 
and dark energy are part of a single component. The key parameter 
in this model is the speed of sound (today) of the dark component, 
c^, that has to be different from zero but small enough to let the 
dark component cluster (see Bertacca et al. 2008, 201 1 for details). 
Detecting a non-zero speed of sound would be an indication of a 
non-ACDM universe. 

Using NVSS we are able to see differences from the ACDM 
case from = 10"'^, while the plot shows that using the com- 
bined full sky EMU-l-WODAN will allow us to constrain values of 
the sound of speed of — 10^'*, using the ISW effect. Further 
details and forecasts on how well we will be able to test UDM cos- 
mological models with these surveys will be part of a subsequent 
paper. 

To predict the significance and constraining power of ISW 
measurements with the forthcoming radio surveys, in Fig. 9 we de- 
fine AC®"^/C®^ as the width of the entire l-cr constraint, i.e.: 

(33) 

where cr^gT is the error on the cross-correlation function in real 
space. In Fig. 9 we show AC"''^/C"'^ (Eq. 33) for the total surveys. 
We compare these with current measurements of AC^^/C^"^ from 
cross-correlations of NVSS and SDSS LRGs with WMAP maps, 
along with the threshold to actually distinguish between ACDM 
and other cosmological models. 

As one can see, for small values of 6 the constraining power is 
maximum, and all the surveys considered should have an increased 
discriminatory power; in the case of whole-sky combined surveys, 
we obtain that AC^^/C^^ is less than half of that of NVSS and 
SDSS, for small 6. 




2 4 6 8 10 12 14 16 18 20 



e (deg) 

Figure 9. Constraining power of cross-correlation CMB-radio sources for 
the different surveys; lines are radio surveys used in this paper, symbols are 
measurements from NVSS and SDSS, lines connecting symbols are thresh- 
olds to detect UDM models with a non-zero speed of sound (see text for 
details). 



5.3 Magnification Bias 

We compute the power spectra for the magnification bias using 
equation (30); we consider the experiment where background ra- 
dio sources from LOFAR, EMU and WODAN are cross-correlated 
with foreground galaxies from SDSS for the northern hemisphere 
and DES for the southern one. We note that we expect to have much 
better and wider data than the SDSS on the timescale of the ra- 
dio surveys considered with the Pan-STARRs Stt sr survey (Kaiser 
et al. 2010), which when complete will provide imaging data to a 
depth between SDSS and DES so our analysis should be consid- 
ered as conservative. In the northern sky we use SDSS galaxies 
up to z = 0.35 and radio sources as the background for higher 
redshift, while for the southern sky we use DES foregrounds up 
\.o z = 1 and EMU sources as a background. To avoid the over- 
lap between foreground and background galaxies, we remove the 
LOFAR, EMU and WODAN galaxies at z < 1, i.e. we assume 
that via cross-matching between optical and radio bands, low z 
radio sources can be removed from our sample. We assume that 
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the bias for the foreground galaxies is unity since they are located 
at low redshifts. We follow Scranton et al. (2005) to measure the 
weighted average power law slope < a — 1 >, where a is given by 
Eq. 18; from our simulations, we obtain -0.219, -0.147, 0.1027 
and 0.121, for LOFAR MS^ LOFAR Tier 1, EMU and WODAN 
respectively. 

In Fig. 10 we show the cross-correlation of radio background 
with optical foreground sources power spectra, computing the er- 
rors according to Zhang & Pen (2006): 



C%g + (Cg + C's'iiot ) (C's + C'shot) 



m + l)/sky 



(34) 



6 COSMOLOGICAL MODEL CONSTRAINTS 

Having calculated the constraints on each cosmological probe, it is 
now possible to determine how the LOFAR, EMU and WODAN 
surveys can improve measurements of cosmological parameters. 
Starting from Einstein's field equations: 



(35) 



where Gp,y is the Einstein tensor, T^ji, is the energy-momentum ten- 
sor and G is Newton's gravitational constant, we study the improve- 
ments these surveys will bring to the measurement of parameters in 
the dynamical dark energy (DE) and modified gravity (MG) sce- 
narios; we investigate these issues using Fisher Matrix techniques, 
following Zhao et al. (2009). 



where / and h denote the foreground and background sources, re- 
spectively, and "shot" stands for the shot noise. 

We can see that using this probe we have better constrain- 
ing power when we have higher number density and magnifica- 
tion index. EMU will provide moderate constraints for this probe in 
combination with DES; we emphasise that a limiting factor for the 
northern surveys (LOFAR and WODAN) in our analysis is actually 
the optical data. We would therefore be able to tighten constraints 
for these surveys if Pan-STARRS data were used as a foreground 
instead of SDSS. 



6.1 Dynamical Dark Energy 

In the cosmological framework of General Relativity, it has been 
necessary to modify Eq. 35 to account for the observed acceler- 
ation of the expansion of the Universe; the simplest modification 
is the introduction of a cosmological constant, as first suggested by 
Zel'dovich (1967), that can be interpreted as vacuum energy; in this 
case one modifies the right hand side of Eq. 35: 

Gf_iv — Tp,y + T™'^, (36) 
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Figure 10. Magnification power spectra for LOFAR-SDSS, EMU-DES and WODAN-SDSS; blaclc solid lines are theoretical predictions (Eq. 30), shaded 
regions are errors, as in Eq. (34). 



where throughout this section we set G = I/Stt and c = 1 for 
simplicity, and: 

= -Ag^,.. (37) 

The first step toward understanding the nature of dark energy is to 
clarify whether it is a simple cosmological constant or it originates 
from other sources that dynamically change in time. The dynami- 
cal models can be distinguished from the cosmological constant by 
considering the evolution of the equation of state of dark energy: 

w=-, (38) 
Q 

where p and g are the pressure density and energy density of the 
fluid, respectively. In the cosmological constant model, w = —1, 
while for dynamical models w = w{a). 

To evaluate the potential of the considered radio surveys to 
constrain the dynamics of different models of dark energy, we adopt 
the following parametrisation for the DE equation-of-state (EoS) w 
(Linder 2003): 

w{a) = wo + Wail — a). (39) 

We use the best fit model {wo^Wa} = {—0.89,-0.24} from 
current data (see Zhao & Zhang (2010) for details) as the fidu- 
cial model, which is consistent with the prediction of the quintom 
model (Feng et al. 2005), and consistently include the dark energy 
perturbations in the calculation using the prescription proposed in 
Zhao et al. (2005). 



6.2 Modified Gravity 

An intriguing alternative to dark energy for the explanation of the 
accelerated expansion of the universe is the "modified gravity" ap- 
proach (Durrer & Maartens 2008), which states that gravity needs 
to be modified, i.e. weakened on large-scales; an attractive feature 
of modified gravity models is that one can alter the Einstein-Hilbert 
action so that accelerated solutions of the background of the Uni- 
verse can be obtained without the need for a dark energy compo- 
nent. 

In this case we modify the geometric side of Eq. 35: 

Gfj,u + G = T^iy , (40) 

Modified gravity models can mimic the ACDM model in the 
sense that they include the background expansion, but in general 
they predict different dynamics for the growth of cosmic structures. 
Radio source number counts and ISW measurements directly probe 
structure formation, therefore they can constrain modified gravity 
scenarios. 

Here we follow Zhao et al. (2010) and consider scalar metric 
perturbations around a FRW background for which the line element 
in the conformal Newtonian gauge is: 

da^ = -a^(r) [(1 + 2*) dr^ - (1 - 2$) df^] , (41) 

where $ and ^ are functions of time and space. 

We use the following parametrisation to describe the relations 
specifying how the metric perturbations relate to each other, and 
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how they are sourced by the perturbations of the energy-momentum 
tensor: 



- = v{a,k), 



-A-KGa^ti{a, k)gA 



(42) 



(43) 



where A is the gauge-invariant comoving density contrast defined 
as: 



A = 5 + 3 — v; 



(44) 



ri{a, k) and ^(a, k) are two time- and scale-dependent functions 
encoding the modifications of gravity that can be written as: 



7?(a,fc) = 



l + Xfk'^a^ 
l + P2\lk^a' 



(45) 



(46) 



where and Pi are parameters and a" gives the time dependence 
of the deviation from GR; rj{a, k) = ^{a, fc) = 1 in GR, while in 
a modified gravity model fi and rj can in general be functions of 
both time and scale (Bertschinger & Zukin 2008; Zhao et al. 2009; 
Linder 2011). 

Since we are interested in testing GR at late times, we will 
consider a simple approximation to Equation 45 and 46 where we 
assume fj,{a, k) = r;(a, fc) = 1 at early times, with a transition 
to some other values at late times. This is natural in the existing 
models of modified gravity that aim to explain the late-time accel- 
eration, where departures from GR occur at around the present day 
horizon scales. Also, the success in explaining the BBN and CMB 
physics relies on GR being valid at high redshifts. 

To model the time evolution of /i and rj we use the hyperbolic 
tangent function to describe the transition from unity to the con- 
stants Ho and rjo: 

K2)=i^(l + tanhi^) (47) 



2 

1 — m 



tanh 



Z — Zs 



+ rjo 



(48) 



2 V Az 
where Zs denotes the threshold redshift where we start to mod- 
ify gravity, and HOjijo are free parameters; following Zhao et al. 
(2010), we fix the transition width Az to be 0.05. 

6.3 Observables and Fisher Matrices 

We use the observables including the LOFAR, EMU and WODAN 
radio source auto-correlation, C^^, cross-correlation, Cf"^ , and 
magnification bias, Cf ^, functions in a Fisher analysis. To obtain 
the auto-correlation functions we consider radio source predicted 
distributions for LOFAR, EMU and WODAN; for the ISW sig- 
nal we cross-correlate the radio source distributions with the CMB, 
while to obtain the magnification correlations we use SDSS DR7 
and DES galaxy populations as foreground lenses for the northern 
and southern hemispheres, respectively. The two-point functions 
we use can be generalised as: 



: An 



Hk 

—A^{k)We''ik)Wr{k), 



(49) 



where A^(fc) is the power spectrum and W^'"^ (k) denote angular 
window functions. Here X,Y £ [T,g, ji], where T, g and fj, indi- 
cate the CMB temperature, radio source counts and magnification 
respectively. 



Given the specifications of the proposed future surveys, the 
Fisher matrix (Fisher 1935, Tegmark et al. 1997) enables us to 
quickly estimate the errors on the cosmological parameters around 
the fiducial values. For Gaussian-distributed observables, such as 



Cf^, the Fisher matrix is given by: 



al3 



/sky 



21+1 



Tr 



dpa ^ dpfj * 



(50) 



where Pa(^) is the a{j3)-lh cosmological parameter and Ct is 



j(/3) 

the "observed" covariance matrix with elements Ci^ 
contributions from noise: 



= Cf ^ + 



that include 



(51) 



Eq. (50) assumes that all fields X(n) are measured over contiguous 
regions covering a fraction /sky of the sky. The value of the lowest 
multipole can be estimated from Imin ~ [7r/(2/sity)], where the 
square brackets denote the rounded integer; for the noise matrix 
N^^' we use Eq. (31, 32, 34). 

To perform the Fisher analysis, we first parametrize our cos- 
mology using: 

P = (tj6,ajc,8a,T,ns, (52) 

where Ub = Q,bh? and Uc = Q,ch^ are the physical baryon and cold 
dark matter densities relative to the critical density respectively, Os 
is the ratio (multiplied by 100) of the sound horizon to the angu- 
lar diameter distance at decoupling, r denotes the optical depth to 
re-ionization, Us and As are the primordial power spectrum index 
and amplitude, respectively, and H £ [woj-qo],!! £ [wa, j^o] are the 
parameters we want to measure. We assume a flat Universe and an 
effective dark energy equation of state w = —1 throughout the ex- 
pansion history; we also combine the latest supernovae la luminos- 
ity distance from the UNI0N2 sample (Amanullah 2010) to tighten 
the constraints. 

Finally, given the uncertainties in the measurement of the bias 
and the redshift distribution for radio surveys, we marginalise over 
the amplitude of the product 6(2) x N{z). We note that the models 
we use are constrained by the total radio source counts and our cur- 
rent knowledge of the evolution of the sub-populations. The main 
uncertainties in these distributions is in the high-redsift (z > 1) 
evolution of the FRI radio galaxies (see e.g. Clewley & Jarvis 2004; 
Sadler et al. 2007), however rapid progress on pinning down the 
evolution of these source should be made over the next few years 
by combining deep multi-wavelength survey data with deep radio 
continuum data (e.g. Smolcic et al. 2009; McAlpine & Jarvis 201 1). 
The final results also depend on the shape of this product, which 
is not precisely known; however, we verified that modifications at 
the level of a few percent in the peak position, amplitude or width 
do not significantly affect our results. A complete analysis of the 
impact of this uncertainty on the measurement of cosmological pa- 
rameters is beyond the scope of this paper and it is left for future 
work. Of course, a careful treatment of this issue will be required 
in the real data analyses, as we mention in Section 8.1. 

We use MGCAMB (Zhao et al. 2009)^ to calculate the observ- 
ables in modified gravity for LOFAR, EMU and WODAN, and 
use Eq. (50) to calculate the Fisher matrices using the preferred 
model of current data as a fiducial model; following Zhao et al. 
(2009) and Zhao et al. (2010), we assume as fiducial {wojWa} = 
{—0.89, —0.24} for the dynamical dark energy parameters, and 
{rjo, fJ-o} = {1-3, 0.87} for the modified gravity parameters. 

^ http://userweb.port.ac.uk/ zhaog/MGCAMB.html 
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7 RESULTS 

The results of our forecasts are shown in Fig. 11, 12, 13, and 14; we 
plot the limits it will be possible to obtain, using the surveys consid- 
ered, in the measurements of the dynamical dark energy and modi- 
fied gravity parameters. To highlight the constraining ability of dif- 
ferent observables, we show the contours for different data combi- 
nations: lighter grey areas are limits from the Planck CMB* plus 
Supemovaela measurements (Amanullah 2010), while darker grey 
areas are improvements we will have adding the auto-correlation 
of radio sources, the radio sources-CMB cross-correlation (ISW), 
the foreground galaxy-background radio source cross-correlation 
(Cosmic Magnification) and a combination of all the measure- 
ments; the crosses refer to the standard model (cosmological con- 
stant and General Relativity), stars indicate the current best fit from 
a combination of probes using WMAP, SDSS and CFHTLS (see 
Zhao et al. 2010 for details). 

We can see that the precision in the measurements of cosmo- 
logical parameters will be significantly increased by the addition of 
the probes considered; in particular we note that the ISW effect is 
more powerful in testing models for gravity than models of dark en- 
ergy. If it turns out that gravity needs to be modified, the ISW effect 
measured with radio surveys will be a powerful probe to measure 
the modified gravity parameters, the physical reason being that if 
/io transits from 1 to another value at low z, indicating a deviation 
from GR, the growth will be enhanced; this will change ($ + 'I') 
significantly, hence generating a large ISW signal (Eq. 13). 

Analyses of clustering and magnification bias also tighten the 
constraints on gravity; the magnification signal measures informa- 
tion about the power spectrum of ($ + ip), which is largely con- 
trolled by 770, and also tests /^o via growth of structures. The ACF is 
also sensitive to ^0 for this reason. These probes will also be useful 
in measuring parameters of the dark energy component, if GR is 
correct even at the largest scales. This is because w{a) changes the 
growth in a very smooth way; so while it does not generate a large 
ISW signal, it does change ($ + 'I') power integrated along the line 
of sight, so can be noticed by magnification bias and the projected 
ACF 

Looking at the different surveys, we can see that they will 
allow precise measurements of cosmological parameters. LOFAR 
Tierl, EMU and WODAN should all be able to increase the preci- 
sion in the dark energy and modified gravity measurements, com- 
pared with that predicted for CMB-l-SNe, by a significant amount; it 
is also interesting to note that adding measurements from LOFAR 
MS'^, which is the least powerful of the surveys we considered (due 
to the lower number density of sources), will already decrease the 
errors in the measurements on modified gravity parameters with 
respect to the CMB-l-SNIa ones. 

Finally, in Fig. 16 we show the constraints on the parameters 
of dynamical dark energy and modified gravity that will be possible 
to obtain using the combination of EMU and WODAN, and we 
compare them to the current best measurements available (Zhao et 
al. 2010). We see that there is substantial improvement, which we 
quantify in Table 2; this reports limits on the measurements of the 
four parameters for the different techniques using the single surveys 
and the EMU-l- WODAN combination. 



http://www.rssd. esa.int/index.php?project=planck 



8 DISCUSSION AND CONCLUSIONS 

In this paper we have presented a forecast of the cosmological mea- 
surements that will be possible with data from the forthcoming LO- 
FAR, EMU and WODAN radio surveys. We have used the correla- 
tion spectra of radio sources: the auto-correlation, the correlations 
with the CMB and with foreground galaxies, alone and in combi- 
nation, to predict measurements of cosmological parameters. 

As mentioned in Section 3.1, the EMU and WODAN surveys 
can be combined in order to obtain a complete full sky catalogue 
and so the largest possible sky coverage, and Fig. 16 shows the 
improvements this combination will produce in constraining cos- 
mological parameters. However, this combination and all the mea- 
surements we highlighted will require a very careful treatment of 
observational data and systematic errors; future work will concen- 
trate on detailed analysis of these issues. 

8.1 Implications for Survey Design 

The tests described in this paper will be very sensitive to system- 
atic errors. For example, to measure magnification bias requires 
that the background source samples are uniformly surveyed (or that 
the threshold variation and completeness are well understood) over 
large areas, placing a stringent requirement on the flux calibration 
of the surveys. Systematics such as these lead to a number of re- 
quirements on the surveys; here we make some initial comments 
about the nature of these requirements: 

• Uniformity and completeness. It is important that the tests 
described in this paper are either conducted on a uniform sample, 
or one where fluctuations are well understood. For instance a 
uniform sample could be created by imposing a flux-density 
cut which is sufficiently above the sensitivity limit at the most 
insensitive part of the survey, so that there are few spurious 
sources, and so that sources are not being lost to systematic effects. 
Detailed simulations will be necessary to check the impact of the 
flux threshold given the consequent non-uniform signal-to-noise. 

• Calibration accuracy of individual surveys. Most surveys 
typically aim for a 1% calibration accuracy; it will be important to 
try to maintain this level, given the need for uniformity described 
above, and the problems arising if these calibration errors occur 
systematically and not randomly across the field. 

• Dynamic Range. If a strong radio source causes low-level 
artefacts, then that will affect the claimed number of faint (and 
therefore typically distant) galaxies, resulting in a spurious cor- 
relation between low-redshift and high-redshift galaxies; to first 
order one will see this as an increase in rms map noise towards 
bright sources. Any of the cosmic measurements need to take this 
into account, possibly through masking the affected area, with the 
consequence of reducing the sky-coverage slightly. 

• Cross-calibration of different surveys. It would be useful for 
all of the surveys to overlap in some regions of the sky to ensure 
an accurate absolute flux scale. 

• Large scale gradients, especially in the declination direction, 
are virtually unavoidable due to changing UV coverage as a 
function of declination and increased system temperatures for low 
elevation observations. These need to be carefully corrected. 

• Bias and redshift distribution uncertainties; this is a well 
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known issue for both the galaxy-galaxy and galaxy-CMB tempera- 
ture spectra and for the redshift distribution only for cosmic magni- 
fication. To take this uncertainty into account, we marginalised over 
the amplitude of b[z) x N{z) (see Section 6.3 for more details). A 
reliable measurement of redshift and bias for the radio continuum 
population will allow us to improve the constraining power of the 
techniques considered in this paper. This is the subject of a future 
paper (Lindsay et al. in prep.). 

8.2 Additional Measurements 

In addition to the techniques presented in this paper, LOFAR, EMU 
and WODAN data will enable several other cosmological analyses, 
which will be useful to test and improve our models. As examples, 
we briefly mention two interesting possibilities: the measurement 
of a dipole anisotropy and a study of the CMB Cold Spot. 

The measurement of a dipole anisotropy in the distribution 
of radio sources can be used to test the distribution of matter at 
different distances and constrain our local motion with respect to 
the comoving cosmic rest frame. 

The dipole anisotropy in the cosmic microwave background 
has been detected with good precision, so an accurate measurement 
of the dipole anisotropy in the large scale mass distribution at lower 
redshift will allow a test of the homogeneity of the matter distribu- 
tion in the universe: if there is agreement between the dipole in 
the CMB and the dipole of galaxies, this will suggest a large scale 
homogeneity; while a discrepancy between the CMB and nearby 
dipole would cast doubt on the general assumption of isotropy and 
homogeneity of the Universe on large scales. 

It is valuable to have radio sky surveys at different frequen- 
cies (such as LOFAR and WODAN), as the amplitude of the radio 
dipole is not only a function of our peculiar velocity, but also of the 
spectral index of radio emission (Ellis & Baldwin 1984). 

A detection of the dipole anisotropy in the radio source distri- 
bution has been reported using NVSS (Blake & Wall 2002), but the 
significance of this measurement depends strongly on the number 
of sources; the surveys considered here will provide an impressive 
improvement in the precision of the dipole anisotropy measure- 
ment, being able to move from an uncertainty of ~15 degrees in 
dipole direction of (Blake & Wall 2002), to an improved accuracy 
of ~2 deg, at 1-cr level (Crawford 2009). 

Using the radio source distribution, it will also be possible 



to perform a number count analysis in order to search for a void 
in the direction of the Cold Spot (Cruz et al. 2005) in the Cos- 
mic Microwave Background. Several models have been proposed 
in order to explain this anomaly, e.g. voids (Inoue & Silk 2006, 
Rudnick 2007), second order gravitational effects (Tomita & Inoue 
2008) or a brane-world model (Cembranos et al. 2008); Cruz et 
al. (2007) showed through a Bayesian statistical analysis that the 
cosmic texture explanation is favoured over the Rees-Sciama effect 
(Rees & Sciama 1968) due to a void or the Sunyaev-Zel'dovich ef- 
fect (Zel'dovich & Sunyaev 1969) caused by a cluster. Radio and 
optical data have been used to test the void hypothesis (Granett et 
al. 2010, Bremer et al. 2010), trying to find a gap in the number den- 
sity in the direction of the Cold Spot; no gap was found, however 
a further analysis using the EMU survey will be helpful, because 
the larger number density of sources at high redshifts will provide 
much better S/N for a potential void. Such an anomalously large 
void will also leave an imprint on ISW measurements (Granett et 
al. 2008), and that again can be examined using EMU data. 

8.3 Conclusions 

In this paper we have shown the potential of SKA pathfinder- 
generation radio surveys to provide competitive cosmological mea- 
surements able to test cosmological models and constrain parame- 
ters describing fundamental physics models. 

Using simulated catalogues, we have predicted which mea- 
surements we will obtain with the source auto-correlation, the 
cross-correlation between sources and the CMB, the magnification 
bias, and a joint analysis together with the CMB power spectrum 
and Supernovae la. 

We have shown examples of the constraining power in testing 
cosmological models alternative to the ACDM-l-GR model, looking 
for modifications coming from non-Gaussianity, alternative models 
for dark energy or modifications to the theory of gravity. We have 
assumed that the surveys will achieve their target dataset and 
treatment of systematic errors, but have tried to be conservative in 
our analyses (e.g. marginalising over the amplitude of correlation 
power spectra, and using objects detected at lOcr signal-to-noise 
threshold). 

There are a number of other galaxy surveys at different wave- 
lengths, which aim to measure cosmological parameters, and which 
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Figure 15. Comparison of median redshift and sky coverage of selected 
future imaging surveys. 



are already collecting data or are being actively prepared for. The 
radio surveys discussed in this paper are complementary to these 
surveys, because of the difference in area, redshift and number den- 
sity covered, and so they will be able to provide useful information 
using some specific probes (i.e. ISW and Cosmic Magnification, as 
their constraining power is increased for larger sky coverage and 
higher redshifts). In the period before SKA, 3D redshift surveys 
such as BOSS will provide more information on the power spec- 
trum on intermediate scales and at low redshifts; photometric sur- 
veys such as Pan-STARRSl^ and DES will also span a different 
part of the parameter space, as they will observe a larger number 
of objects, but at a lower median redshift and, in some cases, a 
smaller region of the sky. Radio surveys cover larger volumes, and 
so provide more large-scale information; thus they will be comple- 
mentary to these other surveys. Next generation experiments such 
as Euclid and LSST will improve the quality of available data, but 
for some aspects the radio surveys of the current generation are 
still competitive, as can be seen from Fig. 15. In the radio, NVSS 
has been used to perform cosmological analyses (e.g. Raccanelli 
et al. 2008, Xia et al. (2010)), and the surveys we considered here 
will have higher median redshift and number of objects observed, 
so they should improve the precision of the cosmological measure- 
ments available. 

Our results show that the unprecedented combination of sky 
coverage, redshift range and sensitivity will enable high-precision 
measurements, competitive with current surveys in a conservative 
scenario. Examining Fig. 16 and Table 2, it is clear that the mea- 
surements that LOFAR, EMU and WODAN could provide are po- 
tentially decisive in ruling out a large part of the cosmological pa- 
rameter space for dark energy and modified gravity models. 
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Table 2. Errors on measurements of dark energy and modified gravity pa- 
rameters for the different surveys and probe combinations (current best 
measurements: {wq, Wa} = {-0.89 ± 0.11, -0.24 ± 0.56} , (Zhao & 
Zhang 2010), {rio, /^o} = {1.3 ± 0.35, 0.87 ± 0.12}), (Zhao et al. 2010). 



of an RCUK fellowship. WJP is grateful for support from the Eu- 
ropean Research Council, the Leverhulme Trust and STFC. DIS 
acknowleges financial support from Deutsche Forschungsgemein- 
schaft (DFG). We thank Daniele Bertacca, Emma Beynon, Annal- 
isa Bonafede, Rob Crittenden, Olivier Dore, Tommaso Giannan- 
tonio, Ben Hoyle, Minh Huynh, Joseph Lazio, Alejo Martinez- 
Sansigre, Sabino Matarrese, Matthias Rubart, Hana Schumacher, 
Charles Shapiro, David Wands and lun-Qing Xia for valuable 
discussions. Numerical computations were carried out on the 
SCIAMA High Performance Compute (HPC) cluster which is sup- 
ported by the ICG, SEPNet and the University of Portsmouth. 



REFERENCES 

Abazajian K. N., et al., 2009, ApJS, 182, 543 
Abdalla F B., & Rawfings S., 2005, MNRAS, 360, 27 
Abdalla F B., Blake C, Rawlings S., 2010, MNRAS, 401, 743 
AmanuUah R., et al., 2010, Astrophys. J., 716, 712 
Bartelmann M., Schneider P, 2001, Phys. Rep., 340, 291 
Bartolo N., Komatsu E., Matarrese S., Riotto A., 2004, Phys. 

Rept., 402, 103 
Baugh, C. M., 2008, PTRSA, 366, 1884, 4381 
Becker R.H., White R.L., Helfand D.J., 1995, AJ, 1 12, 407 
Bennett C. L., et al., 201 1, ApJS, 192, 17B 
Bertacca D., Raccanelli A., Piattella O., Pietrobon D., Bartolo N., 

Matarrese S., Giannantonio T., 201 1, JCAP, 03, 039 
Bertacca D., Bartolo N., Diaferio A., Matarrese S., 2008, JCAP, 

10, 023 

Bertschinger E., Zukin R, 2008, Phys. Rev. D, 78, 024015 
Bessel F, 1824, Abh. d. K. Akad. Wiss. Berlin, 1 
Blake C, Wall J. 2002, Nature, 416, 150 
Boughn S., Crittenden R., 2002, PhRvD, 88, 021302 
Boughn S., Crittenden R., 2004, Nature, 427, 45 



18 A. Raccanelli et al. 




Bremer M. N., Silk J., Davies L. J. M., Lehnert M. D., 2010, MN- 
RAS, 404, L69 

Cabre A., Fosalba P., Gaztanaga E., Manera M., 2007, MNRAS, 
381, 1347 

Cembranos J. A. R., de la Cruz-Dombriz A., Dobado A., Maroto, 

A. L., 2008, arXiv:0803.0694 

Clewley L., & Jarvis, M. J., 2004, MNRAS, 352, 909 
Cole S., Kaiser N., 1989, MNRAS, 237, 1127 
Condon J., et al., 1998, AJ, 115, 1693 

Copi C. J., Huterer D., Schwarz D. J., Starkman G. D., 2010, 
AdAst, 2010, 83C 

Crawford R, 2009, Astrophys. J., 692, 887 

Crittenden R. G. and Turok N., 1996, Phys. Rev. Lett. 76, 575 

Cruz M., Martinez-Gonzalez E., Vielva P., Cayon L., 2005, MN- 
RAS, 356, 29 

Cruz M., Turok N., Vielva P., Martinez-Gonzalez E., Hobson M., 

2007, Science, 318, 1612 
Dalai N., Dore O., Huterer D., Shirokov A., 2008, PhRvD, 77, 

123514 

De Felice A., Koyama K., Raccanelli A., Tsujikawa S., Zhao G.- 

B. , in preparation 

Desjacques V. and Seljak U., 2010, Class. Quant. Grav., 27, 
124011 

Dupe F. X., Rassat A., Starck J. L., Fadili M. J., 2010, 
arXiv:10102192vl 

Durrer R., Maartens R., 2008, Gen. Rel. Grav., 40, 301 

Eisenstein D., et al., 201 1, arXiv: 1 101.1529 

Ellis G. R R., Baldwin J. E., 1984, MNRAS 206, 377 

Feng B., Wang X., Zhang X., Phys. Lett. B, 2005, 607, 35 

Risher R. A., 1935, J. Roy. Stat. Soc, 98, 39 

Giannantonio T., Crittenden R. G., Nichol R. C, Scranton R., 
Richards G. T., Myers A. D., Brunner R. J., Gray A. G., Con- 
nolly A. J., Schneider D. P, 2006, PhRvD, 74, 063520 

Giannantonio T., Scranton R., Crittenden R.G., Nichol R.C., 
Boughn S.P, Myers A.D., Richards, G.T., 2008, PhRvD, 77, 12 

Giannantonio T, Song Y.-S., Koyama K., 2008, PhRvD, 78, 
044017 

Granett B. R., Szapudi L, Neyrinck M. C, 2010, ApJ, 714, 825 
Granett B.R., Neyrinck M.C., Szapudi I., 2008, ApJ, 683, L99 
Hamilton, A.J.S. 1993, ApJ, 417, 19 

Hildebrandt H., van Waerbeke L., Erben T, 2009, A&A, 507, 683 
Ho S., Hirata C, Padmanabhan N., Seljak U., Bahcall N., 2008, 



PhRD 78, 043519 
Inoue K. T, Silk J., 2006, ApJ, 648, 23 

Johnston, S., et al. 2008, Experimental Astronomy, 22, pp 151 
Jonas J., 2009, Panoramic Radio Astronomy: Wide-field 1-2 GHz 

Research on Galaxy Evolution, 
Kaiser N., et al., 2010, Proc. SPIE, 7733 
Kazin E., et al., 2010, ApJ, 710, 1444 

Komatsu E., 2010, Classical and Quantum Gravity, 27, 124010 
Komatsu E. et al, 2011, Astrophys. J. Suppl., 192, 18 
Landy S.D., Szalay A.S., 1993, ApJ, 412, 64 
Laureijs R., et al., 2009, arXiv:0912.0914 
Legendre A.M., 1785, Mm. Math. Phys., 10, 411 
Linder E.V., Phys. Rev. Lett., 2003, 90, 091301 
LinderE.V., 2011, arXiv:1103.0282vl 
Lindsay S., et al., in preparation 

Loverde M., Hui L., Gaztanaga E., 2007, PhRvD, 043519 
Massardi M., Bonaldi A., Negrello M., Ricciardi S., Raccanelli 

A., De Zotti G., 2010, MNRAS, 404, 532 
Matarrese S., Coles P., Lucchin R, Moscardini L., 1997, MNRAS, 

286, 115 

Matarrese S., Verde L., Jimenez R., 2000, ApJ, 541, 10 
McAlpine K., & Jarvis, M. J., 2011, MNRAS, 413, 1054 
Menard B., Scranton R., Rukugita M., Richards G., 2010, MN- 
RAS, 405, 1025 
Mo H.J., White S.D.M., 1996, MNRAS, 282, 347 
Morganti R., et al., 2010, arXiv: 1001.2384vl 
Moscardini L., Coles R, Lucchin R, Matarrese S., 1998, MNRAS, 
299, 95 

Norris R. P., et al. 201 1, PASA accepted, arXiv:l 106.3219 
Oosterloo T., Verheijen M., van Cappellen W., 2010, The latest on 

Apertif. ISKAP2010 Science Meeting 
Overzier R. A., Rottgering H. J. A., Rengelink R. B., Wilman 

R. J., 2003, A&A, 405, 53 
Papal P, Szapudi I., 2010, arXiv: 1009.0754 
Peebles P.J.E., 1980, The large-scale structure of the universe 
Pietrobon D., Balbi A., Marinucci D., 2006, PhRvD, 74, 043524 
Pogosian L., Corasaniti P.-S., Stephan-Otto C, Crittenden R., 

Nichol, R., 2005, PhRvD, 72, 103519 
Press W., Schechter P, 1974, ApJ, 187, 425 
Raccanelli A., Bonaldi A., Negrello M., Matarrese S., Tormen G., 

de Zotti G., 2008, MNRAS, 386, 2161 
Rassat A., 2009, arXiv:0902.1759 



© 0000 RAS, MNRAS 000, 000-000 



Cosmological Measurements with Radio Surveys 19 



Rees M.J., Sciama D.W., 1968, Nature, 217, 511 

Ross A., et al., 2011, arXiv: 1105.2320 

Rottgering H.J.A., 2003, New Astron.Rev., 47, 405 

Rottgering H.J.A., 2010, PoS(ISKAF2010), 050 

Rottgering H., et al., 2011, arXiv: 1107. 1606 

Rudnick L., Brown S., Williams L., 2007, ApJ, 671, 1, 40 

Sachs R. K., Wolfe A. M., 1967, ApJ, 147, 73 

Sadler E. M., et al., 2007, MNRAS, 381, 211 

Samushia L., Percival W. J., Raccanelli A., 2011, 

arXiv:1102.1014vl 
Schaefer B. M., Douspis M., Aghanim N., 2009, arXiv;0903.4288 
Schlegel D., et al., 2011, arXiv: 1106. 1706 
Scranton R., et al., 2005, ApJ, 633, 589 
Sheth R. K., Tormen G., 1999, MNRAS, 308, 119 
Slosar A., Hirata C, Seljak U., Ho S., Padmanabhan N., 2008, 

JCAP08, 031 
Smolcic v., et al., 2009, ApJ, 696, 24 

Thomas S. A., Abdalla F. B., Lahav O., 2010, arXiv: 1012.2272 
Tegmark M., Taylor A. N., Heavens A. P., 1997, ApJ, 480, 22 
Tomita K., Inoue K.T, 2008, PhRvD, 77, 103522 
Van Waerbeke L., 2009, arXiv:0906.1583 
Wang L., et al., 2011, MNRAS, 414, 596 

Wilman R.J., Rottgering H.J. A., Overzier R. A., Jarvis M.J., 2003, 

MNRAS, 339, 695 
Wilman R. J., Miller L., Jarvis M. J., Mauch T, Levrier R, Abdalla 

F. B., Rawlings S., Klckner H.-R., Obreschkow D., Olteanu D., 

Young S., 2008, MNRAS, 388, 1335 
Wilman R. J., Jarvis M. J., Mauch T, Rawlings S., Hickey S., 

2010, MNRAS, 405, 447 
Xia J.-Q., Viel M., Baccigalupi C, Matarrese S., 2009, JCAP, 09, 

003 

Xia J.-Q., Viel M., Baccigalupi C, De Zotti G., Matarrese S., 

Verde L., 2010, ApJ Letters, 717, L17 
Xia J.-Q., Bonaldi A., Baccigalupi C., De Zotti G., Matarrese S., 

Verde L., Viel M., 2010, JCAP, 1008, 013 
Zel'dovich Ya. B., 1967, Zh. Eksp. & Teor. Fiz. Pis'ma 6, 883 
Zel'dovich Ya. B., Sunyaev R., 1969, Astrophys. Space Sci. 4, 301 
Zhang P and Pen U. L., 2006, MNRAS, 367, 169 
Zhao G.-B., Xia J.Q., Li M., Feng B., Zhang X., 2005, Phys. Rev. 

D, 2005, 72, 123515 
Zhao G.-B., Pogosian L., Silvestri A., Zylberberg J., 2009, 

PhRvD, 79, 8, 083513 
Zhao G.-B., Zhang X., 2010, Phys. Rev. D, 81, 043518 
Zhao G.-B., Giannantonio T, Pogosian L., Silvestri A., Bacon 

D. J., Koyama K., Nichol R. C, Song Y.-S., 2010, PhRvD, 81, 

103510 



© 0000 RAS, MNRAS 000, 000-000 



